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Abstract 



,We study two reciprocal thermal effects in the ferromagnetic semiconductor (Ga,Mn)As by scattering theory: domain 
^ Wall motion induced by a temperature gradient as well as heat currents pumped by a moving domain wall. The effective 
C$ put-of-plane thermal spin transfer torque parameter Pq(5q, which governs the coupling between heat currents and a 
magnetic texture, is found to be of the order of unity. Unpinned domain walls are predicted to move at speed 10 m/s in 
1 temperature gradients of the order 10 K//im. The cooling power of a moving domain wall only compensates the heating 
due to friction losses at ultra-low domain wall velocities of about 0.07 m/s. The Seebeck coefficient is found to be of the 
order 100-500 /xV/K at T=10 K, in good agreement with recent experiment. 
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A magnetic domain wall is a region with a gradual reori- 
entation of the local magnetic moments between magnetic 
domains with different magnetization directions (Fig. [T|). 
The position of a domain wall in a magnetic wire can be 
manipulated by an external magnetic field or an electric 
'current [l], • The latter mechanism makes electronic shift 
registers possible in which the magnetic domains separate 
different bits which are collectively moved by a current 
( "racetrack memory" ) Q . 
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Figure 1: (Color online) A ferromagnetic wire containing a domain 
wall in which the magnetic texture rotates in the transverse xz-plane 
(Bloch wall) [||. 

When a current, which in a metallic ferromagnet is spin 
polarized, traverses a magnetic domain wall, it exerts a 
torque on the domain wall, which can be decomposed into 
an in-plane and an out-of-plane component [lj] . In the adi- 
abatic limit, the in-plane torque component can be under- 
stood easily in terms of conservation of angular momen- 
tum. The out-of-plane component is parameterized by a 



so-called /3-factor [5]. Together with the Gilbert damp- 
ing parameter a, which describes magnetic friction pro- 
cesses, it determines the current-driven domain wall mo- 
bility: an increasing /3-factor increases the mobility, while 
it decreases with increasing a [lj]. 

The /3-factor associated with a voltage gradient is de- 
noted /3 C . In transition metal ferromagnets, /3 C and a are 
found to be of the same order of magnitude (~ 10~ 3 — 
10~ 2 ) @, [l|. These parameters have recently also been 
studied in ferromagnetic semiconductors @, H, la] > in which 
the strong spin-orbit coupling in the valence band is found 
to be responsible for a drastically larger /3 c /a ~ 100 @,[§]. 

In the last few years a new research field called Spin 
Caloritronics, loosely speaking the study of thermal effects 
in spintronics, has emerged leading to renewed interest 
in domain wall motion induced by temperature gradients. 
The effect was first studied a couple of decades ago 10 1, 
and has recently enjoyed a renaissance [ll, 12, 13, 14, 15 1. 
Since currents can be induced by both voltage and tem- 
perature gradients, there are two distinct contributions to 
/3. (1) f3 c linked to a voltage gradient, and (2) the recently 
predicted thermal /3q associated with a temperature gradi- 
ent [HI, [3, [lij] ■ Analogous to electrically induced domain 
wall motion, the ability to move the domain wall by a tem- 
perature gradient is controlled by the /3q / a-ratio IH 14 1 . 
It can arise from ballistic non-adiabaticity in narrow do- 
main walls [l6j ] . Here we discuss the adiabatic and dissipa- 
tive contribution to /3q by microscopic scattering theory 
based on the Luttinger Hamiltonian of the GaAs valence 
band. 

In (Ga, Mn)As, strong spin-orbit interaction compli- 
cates the definitions of spin-dependent conductivites and 
Seebeck coefficients. Instead, we calculate the well-defined 
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quantities PqPq and P C /3 C (Pc and Pq denned below). The 
PqPq {P c Pc ) parameter we calculate, describes the ef- 
fective out-of-plane thermal (electric) spin transfer torque 
parameter experienced by a rigidly moving domain wall, 
i.e. a global parameter depending on the magnetization 
profile, in particular the domain wall length. 

(Ga,Mn)As has a large thermopower with Secbeck coef- 
ficients of the order S ~ 100-300 fj,V/K [17]. Here we cal- 
culate the thermopowers of different alloys and find good 
agreement with experiments, implying a strong particle- 
hole asymmetry in the system. Combined with the large 
P c Pc values found before, we also expect a high PqPq in 
these materials. Here we present the first study of domain 
wall motion induced by a temperature gradient in the fer- 
romagnetic semiconductor (Ga,Mn)As. We find a large 
thermally induced domain wall response with PqPq ~ 1, 
which is the same order of magnitude as P C (3 C - This leads 
to the prediction that domain walls move at a speed of 
10 m / s with temperature gradients as low as 100 K / //m, 
much smaller than what has been predicted for transition 
metals [H |. 

The paper is organized as follows: We first review the 
derivation of Pq(3 c and Pr.0p in terms of the scattering 
matrix following Refs. [1, ITjj . Then PqPq is calculated 
numerically for disordered (Ga,Mn)As, which enables us to 
predict the temperature gradients needed to move domain 
walls. We also discuss the heat pumping by domain wall 
systems H, 



Magnetization dynamics is well described phenomeno- 
logically by the generalized LLG equation [l[ 



m = —7m x H e ff + am x m + t, 



(1) 



where m = M(r, t)/M s is the unit direction vector of the 
magnetization M, H off = -SF[M(r,t)]/SM is the effec- 
tive field given by the free energy functional F[M(r, t)], 
7 is (minus) the gyromagnetic ratio, and a the Gilbert 
damping constant. In conducting ferromagnets, the itin- 
erant quasiparticles are strongly coupled to the magnetiza- 
tion via the exchange interaction. Through this coupling a 
spin current can induce a torque r on a magnetic texture. 
The torque can be generated by a temperature gradient 
AT = T L - T R (t q ) or a voltage gradient AV = V L - Vr 
(t c ). To lowest order in magnetization gradients these two 
torques have the form [13J, [l4j 



T c = -f^-AV(l - j3 c mx) dym, 



HPnTSG AT ,„ 
Tn = =r- 1 



(2) 
(3) 



where So — A/s/7 is the spin density along — m, G is 
the conductance, S = —eCTdE^G is the Seebeck co- 
efficient, C = (fcs7r) 2 /(3e 2 ) is the Lorenz constant and 
8e = d/dE at the Fermi energy, e the modulus of the 
electron charge, 2P C = (a^ — a^) / (a^ + a^) where crlw are 
the spin-dependent conductivities, Pq = P c +P s (l — 4P 2 ) 
lj| with 2P S = (S 1 " - S l )/(S^ + S±) where are the 



spin-dependent Seebeck coefficients, and A is the cross sec- 
tion area of the conductor. The first terms in Eqs. © and 
([3]) are the in-plane torques, which describe the angular 
momentum transfer on m by a spin current following the 
magnetic texture adiabatically. The second terms, param- 
eterized by P c and /3q, describe dissipative out-of-plane 
torques. /3q, f3 c and a govern the domain wall response to 
an applied temperature or voltage gradient. 

We adopt a magnetic free energy functional for a wire 
along the y-axis 

F[M] = M s Jdv(^[(V9) 2 + sin 2 (0)(V0) 2 ] + 

^ sin 2 9 sin 2 <$> - ^ cos 2 9 - H cxt cos ^ , (4) 

where K s is the spin wave stiffness, K± and K z are 
anisotropy constants, H cxt is an external magnetic field, 
and 9 and <f> are the polar and azimuth angles describ- 
ing the local magnetization directions, respectively. A 
local minimum of this functional is a Bloch wall [J] ro- 
tating in the (transverse) x—z plane with local magne- 
tization direction given by cos# = tanh([y — r w ]/X w ), 
sin# = 1/ cosh([y — r w ]/X w ), where r w is the position of 
the wall and A^, the domain wall length. The position can 
be manipulated either by a magnetic field or a spin cur- 
rent. For low current densities and magnetic fields the wall 
moves rigidly with constant tilt angle <f> and domain wall 

length \ w = K S /(K Z + Kj_ sin 2 (</>)), while the domain 
wall position r w is given by: 
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Eq. ([5]) is obtained by substituting the Bloch wall ansatz 
with a time-dependent tilt angle <f>(t) into Eq. (fTJ), and 
consider the regime where d(j)/dt = 0. The values of 
AT/T, AV and H cxt at which this regime breaks down 
are known as the Walker thresholds, above which also the 
parameters A„, and 4> have to be treated as dynamic vari- 
ables. In this paper we will always assume AT/T, AV and 
ff cx t to be below the threshold fields, so that r w (t) is the 
only time-varying parameter. In (Ga,Mn)As the Walker 
ansatz breaks down for velocities typically of the order 10 
m/s 0. 

Eq. (J5J) shows that the magnetic system responds to tem- 
perature and voltage gradients. From Onsager's reciprocal 
relations [l9[ we know that the reciprocal processes heat 
and charge pumping by a moving domain wall also exist, 
and that the response coefficients have to obey Onsager 
symmetry relations. The response coefficients for heat 
and charge pumped by a moving domain wall thus gives 
us expressions for the parameters PqPq and P C /3 C - This 
strategy is used in Refs. to derive scattering matrix 

expressions for P C /3 C and PqPq- The response coefficients 
Lq w (L cw ) and L w q (L wc ), describing heat (charge) cur- 
rent pumped by a moving domain wall and domain wall 
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motion induced by a temperature (voltage) gradient, re- 
spectively, are found from the parametric pumping for- 
mula [If SEl: 
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Here, S is the scattering matrix, Tr the trace over prop- 
agating modes in the leads, and f 2 = [l 0; — l] a 
2N x 2N block-diagonal matrix, where 1 is the N x N 
unit matrix and N is number of propagating modes at 
the Fermi energy in each lead, and we assume a mirror 
symmetric system. Eq. © and Q also assume a suffi- 
cient long wire where back-action of the induced currents 
on the domain wall motion can be neglected 14|. On- 
sager's reciprocal theorem implies the symmetry relations 
L Qw [m] = L wQ [-m] and L cw [m] = L wc [-ra]. Eq. © 
and ([7]) give the following expressions for P C (3 C and Pq/3q: 
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(8) 
(9) 



where t is the transmission matrix. 

In addition to the heat current pumped between the 
left and right contacts by the moving domain wall, heat 
is also generated by the magnetization damping. For slow 
magnetization dynamics, this energy dissipation rate can 
be obtained from the scattering matrix as well [2l], HJ : 
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dr w dr n 



(10) 



The damping is second-order in the magnetization dynam- 
ics and therefore not controlled by Onsager's theory. 

We neither took into account inelastic electron scatter- 
ing nor temperature dependence of the material param- 
eters. This is believed to be a good approximation for 
the highly conducting transition metals, but might become 
significant in magnetic semicondutors, but a more detailed 
study is beyond the scope of this manuscript. 

In the following, we focus on the (Ga,Mn)As system, 
and calculate Pq/3q using Eq. © and P C /3 C by the methods 
reported in Ref. |8[- We investigate the reciprocal effects: 
Heat pumping by a moving domain wall as well as torque 
induced by a temperature gradient. 

To model the band structure of (Ga,Mn)As we use the 
Hamiltonian [22J 



H = H L + h(r)-J + V(r), 



(11) 



where Hl is the 4x4 Luttinger Hamiltonian for zincblende 
semiconductors in the spherical approximation. J is a vec- 
tor of 4x4 spin matrices for J— 3/2 spins. The h- J term is 
a mean field approximation of the exchange interaction be- 
tween the itinerant holes and the local magnetic moment 
of the Mn dopants. The exchange field h is antiparallel 
to the magnetization direction m. V(r) = J2i — R-i) 
is the impurity scattering potential, where Rj is the posi- 
tion of impurity i, and Vi its scattering strength (23[. The 
Vi are randomly and uniformly distributed in the interval 
[-V Q /2,V /2]/ 

The response coefficients L cw and Lq w are calculated 
from the scattering matrix expressions in Eq. ^ and 
Eq. ([7])- The scattering matrix is calculated numeri- 
cally using a stable transfer matrix method j24|. Disor- 
der effects are included fully and nonperturbatively by 
the ensemble average = Yln=i (^v)n /^ii where 

Ni is number of different impurity configurations. All 
coefficients are averaged until an uncertainty 5 = 



of less than 10% is achieved. 



The mean free path I for impurity strength Vq is cal- 
culated by fitting the average transmission probability 
T = (G)/G sh to T(Ly) = 1/(1 + L v ) [H, where G sh 
is the Sharvin conductance and (G) is the conductance for 
a system of length L y . 

We consider a discrete (Ga,Mn)As system with trans- 
verse dimensions L x — 19 nm, L z = 15 nm and L y € 
[100, 350] nm connected to infinite ballistic leads. The lat- 
tice constant is 1 nm, much less than the Fermi wavelength 
Xp ~ 10 nm. The Fermi energy is 0.077 eV when measured 
from the lowest subband edge. The Luttinger parameters 
in H L are 7i = 7.0 and 72 = 2.5, and |h| = 0.032 eV j^. 
X w G {10,20,40} nm are the wall widths. To estimate 
a typical saturation value of the magnetization we use 



M s = 10|7|foc/aQ aAs 26], with x — 0.05 as the doping 



level, and where acaAs is the lattice constant for GaAs. 

Since we used the Sommerfeld approximation, the scat- 
tering matrix expressions for the Onsager coefficients in 
Eq. © and Eq. (JTJ) are valid only for thermal energies 
that are small compared to the Fermi energy. For instance 
ksT/EF = 0.01, implies a temperature of the order 10 K 
with the Fermi energy specified above. 

Fig. [5] (a) and (b) show the out-of-plane torque param- 
eters P c fic and PqPq as a function of disorder. We ob- 
serve that in the diffuse regime (I <^ X w ) PqPq ~ 1.0 is 
of the same order as P c f3 c . As explained in Ref. @, the 
high P c f3 c of (Ga,Mn)As arises from large hole reflection 
at the domain wall due to the intrinsic domain wall re- 
sistance caused by spin-orbit coupling p7| . This intrinsic 
resistance decreases with disorder, reducing P c j3 c with in- 
creasing disorder. For even smaller mean-free paths, P C (3 C 
saturates or increases slightly due to the increasing spin- 
flip rate caused by the impurity scattering. In the ballistic 
limit (I ^> X w ) Pq/3q is nearly one order of magnitude 
larger than P C /3 C , but decreases rapidly when disorder is 
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Figure 2: (a) P C /3 C and (b) PqPq plotted as a function of domain 
wall length (X w ) over mean free path (/). The lines are guides to the 
eye. 



included, saturating in the diffusive regime. Since PqPq is 
a measure of the effect of particle-hole asymmetry on the 
out-of-plane torque efficiency, this means that the particle- 
hole asymmetry of the torque efficiency is stronger in the 
ballistic than in the diffuse regime. 

In the diffuse regime we find Sccbcck coefficients S ~ 
fOO - 500 a*V/K (S/T= 10 -50 J^V/K 2 ) at T = 10 K 
for mean free path I ~ 7 — 27 nm [281 ] . in good agreement 
with the experimental values S ~ 100 — 300 /j,V /K 17 1. 
The Gilbert damping calculated from Eq. (|10l) is of the 
order 10" 3 - 10" 2 If. 

The high Pq(3q combined with a large Seebeck coeffi- 
cient implies a strong coupling between thermally induced 
currents and domain wall motion in (Ga,Mn)As. This in- 
creases the ability to manipulate the domain wall position 
with a temperature gradient, and oppositely, to transfer 
heat between the left and right contacts by moving the 
domain wall. These two effects are illustrated in Fig. [3J 
where we plot the reciprocal thermal effects: Domain wall 
motion induced by a temperature gradient and heat cur- 
rent pumped by a moving domain wall. Fig. [3] (a) shows 
the domain wall velocity divided by T as a function of 



temperature gradient, while Fig. [3] (b) shows the "ther- 
mal motive force" Vq — Jq /G, a system-size independent 
quantity determining how efficient the domain wall system 
pumps a heat current in the absence of a temperature or 
voltage bias, divided by T 2 as a function of domain wall 
velocity 
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Figure 3: (a) Plot of domain wall velocity divided by temperature in- 
duced by the temperature gradient (T^ —Tn)/L y . (b) The "thermal 
motive force" Vq = Jq /G is displayed as a function of domain wall 
velocity. In both plots the Seebeck coefficient is of the order S/T = 
10 - 50 jUV/ K 2 and the conductivity a = 270 - 630 Q' 1 cm' 1 . 

To get a better feeling of the importance of these two 
reciprocal processes, let us consider a wire at tempera- 
ture T = 10 K with L y = 60 /^m, and conductivity a ~ 
270 fP 1 cm -1 , containing a Bloch wall with \ w = 10 nm. 
We limit attention to the regime in which the domain wall 
moves rigidly, i.e. below the Walker thresholds. We see 
from Fig. [3] (a) that a temperature gradient of around 
10 K / /xm is needed to drive the domain wall to speed 10 
m/s. According to Eq. ([S]), an electric field-induced cur- 
rent density of the order 10 4 A / cm 2 will give an equiv- 
alent torque. However, in a real (Ga,Mn)As system, we 
know that current densities of the order 10 5 — 10 6 A / cm 2 
are needed to achieve such an effect [29j]. The disagree- 
ment is caused by reduction of the domain wall mobility 
by extrinsic pinning effects [30j. Thus, taking into account 
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this effect, one may expects that a temperature gradient 
of the order 100 — 1000 K / /xm needed for driving the wall 
at speed 10 m / s is a more realistic estimate. 

As mentioned above, there is also Joule heat genera- 
tion associated with the magnetization damping. Note 
that this heating does not lead to a temperature difference 
between the left and right contact for a mirror symmet- 
ric system since an equal amount of energy is pumped 
into the right and left contact, but equally heats both 
two contacts. Assuming a Gilbert damping of the or- 
der a — 10~ 2 implies that the domain wall dissipates a 
heat current density of 30 W/m 2 due to magnetic fric- 
tion Q into each reservoir, while the moving domain wall 
pumpes 0.2 W/m 2 between the left and right contact. 
Thus, close to the Walker threshhold the effect of the 
magnetic heat pump is rather small compared to the mag- 
netic friction process. The heat current density pumped 
by the moving domain wall (j'q) scales linearly with f w , 
while the energy current density associated with magneti- 
zation damping (jg as ) scales quadratically. For instance, 
for the system considered above j'q = 0.02^ sW/m 3 and 
jq as = 0.3^ s 2 W/m 4 . The domain wall system is there- 
fore not an efficient cooler except for very small domain 
wall velocities (r w < 0.07 m/s), where the domain wall 
pumped energy will dominate over the magnetic friction 
process. The situation is the same in transition metal fer- 
romagnets, in which the out-of-plane torque parameter is 
of the same order as the Gilbert damping [1]. Insulating 
ferromagnets such as the Yttrium-Iron-Garnets have very 
low Gilbert damping and are possibly better suited for 
cooling purposes. 

In conclusion, we have studied thermally induced do- 
main wall motion , and the reciprocal effect, a heat pump 
operated by moving the domain wall, in the ferromag- 
netic semiconductor (Ga,Mn)As. The PqPq parameter 
that governs both effects is found to be of the order unity 
in diffuse systems. We estimate that a domain wall veloc- 
ity of 10 m/s can be induced by a temperature gradient of 
the order 100 — 1000 K / fim, while a domain wall pumps 
around 0.2 W / m 2 when it moves at speed 10 m/s. The 
energy dragged by the domain wall from one contact to 
the other is relatively small compared to energy dissipa- 
tion associated with magnetic friction at these velocities, 
implying that (Ga,Mn)As domain wall systems are not 
useful for cooling purposes. 
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